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ABSTRACT 

Aims. We obtained phase-resolved spectroscopy of the accreting millisecond X-ray pulsar SAX J1808.4-3658 during its outburst in 
2008 to find a signature of the donor star, constrain its radial velocity semi-amplitude (A'2), and derive estimates on the pulsar mass. 
Methods. Using Doppler images of the Bowen region we find a significant (>8(t) compact spot at a position where the donor star is 
expected. If this is a signature of the donor star, we measure A"j,„=248±20 km s"' (Icr confidence) which represents a strict lower limit 
to A'2. Also, the Doppler map of He II /14686 shows the characteristic signature of the accretion disk, and there is a hint of enhanced 
emission that may be a result of tidal distortions in the accretion disk that are expected in very low mass ratio interacting binaries. 
Results. The lower-limit on K2 leads to a lower-limit on the mass function of /(Mi)>0.10Mo. Applying the maximum A'-correction 
gives 228<A:2<322 km s"' and a mass ratio of 0.051<^<0.072. 

Conclusions. Despite the limited S/N of the data we were able to detect a signature of the donor star in SAX J1808. 4-3658, although 
future observations during a new outburst are still warranted to confirm this. If the derived Ke,,, is correct, the largest uncertainty in the 
determination of the mass of the neutron star in SAX JI808.4-3658 using dynamical studies lies with the poorly known inclination. 

Key words, accretion, accretion disks - stars:individual (SAX J1808.4-3658) - X-rays:binaries. 



1. Introduction 

Low-Mass X-ray Binaries (LMXBs) are systems in which a 
compact object (a neutron star or a black hole) is accreting mat- 
ter, via Roche lobe overflow, from a low mass (<lMo) star. Some 
. LMXBs exhibit sporadic outburst activity but for most of their 
' time remain in a state of low-level activity (White et al. 1984); 
we will refer to these systems as transients. In April 1998, a co- 
herent 2.49 ms X-ray pulsation was discovered with the Rossi 
X-ray Timing Explorer (RXTE) satellite in the transient SAX 
J1808.4-3658 (Wijnands & van der Klis 1998, Chakrabarty 
& Morgan 1998). This was the first detection of an Accreting 
Millisecond X-ray Pulsar (AMXP). Seven more of these sys- 
tems have been discovered since then; all these systems are tran- 
sients, have orbital periods in the range between 40 min and 4.3 
hr and spin frequencies from 1.7 to 5.4 ms. These findings di- 
rectly confirmed evolutionary models that link the neutron stars 
of LMXBs to those of millisecond radio pulsars via the spinning 
up of the neutron star due to accretion during their LMXB phase 
(e.g. Wijnands 2006). 

To date, six episodes of activity were detected from SAX 
J 1808.4-3658 with a 2-3 year recurrence cycle. During the 
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1998 outburst a detailed analysis of the coherent timing be- 
haviour showed that the neutron star was in a tight binary sys- 
tem with a 2.01 hr orbital period (Chakrabarty & Morgan 1998; 
Hartman et al. 2008). The mass function derived from X-ray 
data (4 x 10"^ Mq) and the requirement that the companion fills 
its Roche lobe led to the conclusion that it must be a rather low 
mass star, possibly a brown dwarf (Chakrabarty & Morgan 1998; 
Bildsten & Chakrabarty 2001). 



We present here optical spectroscopy of SAX J1808. 4-3658 
obtained during the 2008 outburst. One of the key issues in dy- 
namical studies is to measure the radial velocity of the compan- 
ion star (^^2) and use this value to constrain the optical mass 
function of the system. During quiescence, these goals can be 
achieved by tracing the absorption features originating in the 
photosphere of the companion star. However, the intrinsic faint- 
ness of the low mass companion stars in AMXPs makes such 
an analysis very difficult. To overcome this problem, Steeghs & 
Casares (2002) have shown that during phases of high mass ac- 
cretion rates, Bowen blend lines emitted by the irradiated face of 
the companion star can be used. A precise measurement of K2 
represents the only way to determine the optical mass function 
of the system and ultimately constrain the neutron star mass. 
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Fig. 1. Average flux calibrated spectrum of SAX J1808.4-3658. We have labefled the most important features that are present. The 
strong absorption features around 3940 and 3980 A are due to difiiise interstellar bands. 



2. Observations and data reduction 

Optical spectroscopic observations of SAX J1808. 4-3658 were 
carried out on 27 September 2008 with the ESO Very large 
Telescope (VLT), using the 1200B grism on FORSl with a slit 
width of 0.7 arcsec. We obtained 16 spectra of 360 seconds inte- 
gration each, which corresponds to one orbital period. The see- 
ing during the observations was in the range 0.8" - 1.1". 

Image reduction was carried out following standard proce- 
dures: subtraction of an averaged bias frame, division by a nor- 
malised flat frame. The extraction of the spectra was performed 
with the ESO-MIDA$] software package. Wavelength and flux 
calibration of the spectra were achieved using a helium-argon 
lamp and by observing spectrophotometric standard stars. Our fi- 
nal reduced spectra have a wavelength range from 3600-5000 A, 
a dispersion of 0.72 A pixel ' and resolution of R=2200. Cross 
correlation of the spectral lines and the Doppler tomograms were 
obtained usingthe MOLLY and DOPPLER packages developed 
by Tom Marslfl 

3. Data analysis 

We present the average spectrum of SAX J1808. 4-3658 in 
Fig.[T] The spectrum is dominated by strong Balmer lines in ab- 
sorption, which are due to the optically thick disk in the high 
state and are a typical signature of a low-to-intermediate in- 
clination system. The He II /14686 and the Bowen complex (at 
/1/14630-4660) are also clearly detected as emission features, and 
we have indicated the most important lines. 

For other bright LMXBs, narrow components in the Bowen 
emission have been reported that are thought to arise on the ir- 
radiated surface of the donor star ( e.g. Steeghs & Casares 2002; 
Casares et al. 2006; Cornelisse et al. 2007; see also Cornelisse 
et al. 2008 for an overview), and here we attempt to find similar 
features in SAX J 1808. 4-3658. To calculate the orbital phase for 
each spectrum we used the recent ephemeris by Hartman et al. 
(2008), but added 0.25 orbital phase to their phase zero so that it 
represents inferior conjunction of the secondary. Unfortunately, 
the individual spectra do not have sufficient S/N to identify the 
narrow features, and we must resort to the technique of Doppler 
tomography (Marsh & Horne 1988). This technique uses all the 
spectra simultaneously to probe the structure of the accretion 
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disk and identify compact emission features from specific loca- 
tions in the binary system. However, for this technique to work, 
an estimate of the systemic velocity of the system, y, is crucial. 
We do want to note that by changing the phase zero of the ac- 
curate Hartman et al. (2008) ephemeris any potential donor star 
feature in the map must now lie along the positive y-axis in the 
Doppler tomogram. Thus finding a significant feature there will 
give strong support to an emission site located on the irradiated 
donor star 

To find y we started by applying the double-Gaussian tech- 
nique of Schneider & Young (1980) to Hell A4686. Since the 
wings of the emission line should trace the inner-accretion disk, 
it should not only give us an estimate of the already known radial 
velocity of the compact object (Chakrabarty & Morgan 1998), 
but also y. Using a Gaussian band pass with FWHM of 400 km 
s"' and separations between 500 and 1800 km s"' in steps of 50 
km s"' , we find that between a separation of 1400 to 1600 km s"' 
our values for Ki are close to the one obtained by Chakrabarty & 
Morgan (1998), while at larger separations we are reaching the 
end of the emission line (see Fig.|2]i. Also in this range, our fit- 
ted orbital phase zero (<po) is close to 0.5, further suggesting that 
we are tracing the radial velocities of a region close to the neu- 
tron star, while the systemic velocity is stable around -50 km s"' . 
Despite the large errors on our fits, we do think this test already 
gives a good first estimate of y around -50 km s"^ 

Another test to obtain y is to create Doppler maps for He II 
/14686. Contrary to the Bowen region (which is very complex 
due to the presence of many different lines). Hell is a single 
line that is usually a good tracer of the accretion disk (see e.g. 
Cornelisse et al. 2007; Casares et al. 2006). We searched for y 
between -160 and km s ' in steps of 20 km s ', and all the 
maps show the expected accretion disk structure (see Fig.|3]l. We 
must unfortunately conclude that He II is not very sensitive to 
y, and only suggests a range between -160 and km s"'. We 
note that the maps are dominated by an emission feature in the 
top-left quarter of the map, which we interpret as the gas stream 
impact point. We also note that further downstream there is en- 
hanced emission, that might be due to matter streaming along the 
edge of the disk as was observed in for example EXO 0748-676 
(Pearson et al. 2006). Finally we note that there is some en- 
hanced emission in the top-right corner which might be due to 
strong tidal interaction (see below). 

Cornelisse et al. (2008) have shown that the strongest narrow 
component in the Bowen emission is usually NIII /14641. Our 
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Fig. 2. The derived fit parameters of the radial velocity curve of 
the Hell /14686 emission as trace by a double Gaussian with 
separation a. The free parameters are the orbital phase zero, (po, 
compared to the ephemeris of Hartman et al. (2008), its radial 
velocity, Ki , and its systemic velocity, y. The dotted line in the 
top panel indicates the expected phase zero for the compact ob- 
ject. 



next step was therefore to create Doppler maps of the Bowen 
region for y between and -120 km s"' (again in steps of 20 km 
s ') including only this line. Only when y was between -80 and 
-20 km s ' was a clear spot present and centred on the jt:=0 axis. 
For this range we estimated the peak value and FWHM of the 
spot as a function of y, and in Fig.|4]we show how the ratio of 
these values change. Around 7=-50±15 km s"' (Icr confidence) 
the FWHM/peak value reaches a minimum suggesting that here 
the spot is most compact, and this is the value we adopt for the 
systemic velocity. Furthermore, we also noted that in the range 
from -35 to -65 km s"' the velocity centroid of the spot was 
stable between yv=240-260 km s"' . 

To decrease the noise present in the Bowen Doppler map 
we included the most important other lines that are most often 
present in other LMXBs (NIII ^4634 and CIII ^4647/4650), 
and present this map in Fig. [3] To estimate the significance of the 
compact spot we measured the standard deviation of the bright- 
ness of the pixels in the background. We find that the central 
pixels of the compact spot are 18cr above the background, and 
even 8cr above the second most prominent feature in the map, 
namely the one at (-150,-200) in the Bowen map of Fig. [3] (for 
which it is unclear if it is real or an artifact of the tomogram). 
This strongly suggest that the compact spot is real. 

Finally, to optimise our estimate for Ki,„, we created average 
spectra in the rest frame of the donor star, changing Kg,,, in steps 
of 2 km s ' within our error range. We find that N III A464-Q is 
most pronounced for Kem=24S±2Q km s~' (Icr confidence), and 
adopt this as our final value, and in Fig.|5] show the final aver- 
age spectrum in the rest-frame of the donor star We do note that 
this value is smaller than the ^300 km s"' obtained by Ramsay 
et al. (2008) from the same dataset, since they provide no er- 
rors we cannot tell if this difference is significant. However, as 
stated above, the location of the spot remains within the quoted 
error range for a range of assumed y velocities. To search for 
variability in the Bowen lines as a function of orbital phase, we 
created an average spectrum in the rest-frame of the donor us- 
ing only spectra taken between orbital phases 0.25 and 0.75, and 
another corresponding to phases 0.75 and 1.25. The strength of 
N III /J4640 did not change between these spectra, but it is un- 
clear if this is real (suggesting that the inclination is low) or due 
to the limited S/N of the dataset. 
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Fig. 3. Doppler maps of He II /14686 (top) and the Bowen com- 
plex (bottom). Indicated on both maps is the Roche lobe, gas 
stream leaving the Li point, and the Keplerian velocity along the 
stream for the assumptions ^=0.059 and K2-215 km s"' . 



4. Discussion 

We have presented phase-resolved spectroscopy of the accret- 
ing millisecond X-ray pulsar SAX J1808. 4-3658, and detected 
a compact feature in the Doppler map of the Bowen complex. 
Although the S/N of the data is limited, this spot is the most 
stable and significant feature for a large range of y velocities. 
Furthermore, thanks to the very accurate ephemeris (Hartman 
et al. 2008), the spot is at a position where the donor star is 
expected, and we conclude that it is real. Therefore, following 
detections of a donor star signature in other X-ray binaries (e.g. 
Steeghs & Casares 2002; Casares et al. 2006; Cornelisse et al. 
2007; see Cornelisse et al. 2008 for an overview), we also iden- 
tify this feature as being produced on the irradiated surface of 
the donor star We note that in Fig. |5] most peaks in the Bowen 
region appear to line up with known NIII and CIII lines (e.g. 
Steeghs & Casares 2002) when using our derived values for y 
and Kem- 

Since the donor star surface must have a lower velocity 
than the centre of mass, the observed Kem^^A^+lO km s"' (Icr 
confidence) is a lower limit on the true K2 velocity. However, 
it still gives us a strict lower limit on the mass function of 
/(M)=Misin^//(l-H^)^>O.lOM0, where q is the binary mass ratio 
M2/M1 and / the inclination of the system. We can further con- 
strain the mass function by applying the so-called /T-correction 
(Munoz-Darias et al. 2005), and using the fact that /ri = 16.32 
km s"' (Chakrabarty & Morgan 1998). The largest /T-correction 
possible is when we assume that there is no accretion disk and 
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Fig. 4. The ratio of the FWHM over the peak value for the com- 
pact spot in the Bowen map (using only N III /14661) as a func- 
tion of the systemic velocity (y). 



almost all radiation is produced in the Li point. Applying the 
polynomials by Munoz-Darias et al. (2005) for A'(.„,=248±20 km 
s ' gives /r2=299±23 km s which should be independent of 
the inclination of the system. 

This gives conservative estimates of 228<^r2<322 km s"' 
and 0.05 1<^<0. 072, which we used to create the Roche lobe 
and gas streams on the Bowen Doppler map in Fig. [3] We do 
note that our obtained mass ratio is rather extreme, and would 
suggest that tidal interaction in SAX J1808. 4-3658 is important 
enough to produce a precessing accretion disk and thereby a su- 
perhump (see e.g. O'Donoghue & Charles 1996). This might be 
the explanation of the enhanced emission in the top-right quarter 
of the Hell map (Fig.[3]'top), which was for example also ob- 
served in LMCX-2 (Cornelisse et al. 2007). Such a superhump 
should be moving through the accretion disk on the precession 
time scale, and therefore changes position in the Doppler map 
over time. Unfortunately, since we only have one orbit of data 
we cannot test this, and future observations will be needed to 
see if the spot is long-lived and moves, in order to confirm the 
presence of the superhump. 

Despite our large range of K2 we nonetheless will review the 
implications for the mass of the neutron star. First of all we can 
improve our estimate of K2 by taking into account the results by 
Meyer & Meyer-Hofmeister (1982). They analysed the effects 
of X-ray heating on the accretion disk and found that there is a 
minimal disc opening angle of ^6°, even in the absence of irradi- 
ation. Using this value to estimate the ^T-correction, the polyno- 
mials by Munoz-Darias (2005) suggest that /r2<3 10 km s ' , still 
comparable to the maximum correction possible. Deloye et al. 
(2008) used photometry of SAX J1808. 4-3658 in quiescence to 
constrain the inclination between 36 and 67 degrees. Using these 
extreme values for their inclination and 228<A'2<322 km s 
leads to a neutron star mass between 0. 15 and 1 .58M0. Although 
these values are not very constraining for the neutron star mass, 
they do favour a mass near the canonical value rather than a mas- 
sive neutron star. We also note that our values are smaller than 
the >1.8M0 estimated by Deloye etal. (2008)(fora 10% error in 
the distance estimate), but since this corresponds to a 1.7cr dif- 
ference (only taking into account our errors), we conclude that 
this disagreement is marginal. 

5. Conclusions 

The observations presented here provide evidence for the 
detection of the signature of the irradiated donor star in 
SAX J1808. 4-3658. Clearly a similar experiment, but at higher 
S/N and spectral resolution, must be carried out again during 
a future outburst to obtain more than a single orbital period of 
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Fig. 5. Blow-up of the Bowen region for the average spectrum 
(top) and the average in the rest-frame of the donor star (bot- 
tom). The bottom spectrum clearly shows the narrow compo- 
nents. Indicated are the most important N III (A4634/4640) and 
C III (/14647) Unes. 



data and resolve the narrow lines. This will not only allow us 
to unambiguously claim the presence of narrow components in 
the spectra of SAX J1808. 4-3658, but also measure the rota- 
tional broadening of the narrow components to further constrain 
K2 via the relation in Wade & Horne (2003). With these data 
we have shown a promising way forward to constrain K2 and in 
combination with more quiescent data, we should be able to bet- 
ter constrain the inclination and thereby obtain the mass of the 
neutron star in SAX J1808.4-3658. 
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